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Chapter 1 

Ultrafast control of strong-field electron 
dynamics in solids 


Vladislav S. Yakovlev^’^, Stanislav Yu. Kruchinin^ Tim Paasch-Colberg*, 
Mark I. Stockman^, Ferenc Krausz*’^ 


Abstract We review theoretical foundations and some recent progress related to 
the quest of controlling the motion of charge carriers with intense laser pulses and 
optical waveforms. The tools and techniques of attosecond science enable detailed 
investigations of a relatively unexplored regime of nondestructive strong-held ef¬ 
fects. Such extremely nonlinear effects may be utilized to steer electron motion 
with precisely controlled optical helds and switch electric currents at a rate that is 
far beyond the capabilities of conventional electronics. 


1.1 Introduction 

It has long been realized that intense few-cycle laser pulses provide unique con¬ 
ditions for exploring extremely nonlinear phenomena in solids [1, 2], the key idea 
being that a sample can withstand a stronger electric held if the duration of the 
interaction is shortened. Ultimately, a single-cycle laser pulse provides the best con¬ 
ditions for studying nonperturbative strong-held effects, especially those where the 
properties of a sample change within a fraction of a laser cycle. The recent rapid 
development of the tools and techniques of attosecond science [3] not only creates 
new opportunities for detailed investigations of ultrafast electron dynamics in solids, 
but it also opens exciting opportunities for controlling electron motion in solids with 
unprecedented speed and accuracy. Conventional nonlinear phenomena that accom¬ 
pany the interaction of intense laser pulses with solids have already found a vast 
number of applications in spectroscopy, imaging, laser technology, transmitting and 
processing information [4]. It can be expected that the less conventional nonpertur- 
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bative nonlinearities may also find important applications [5, 6]. The main purpose 
of this chapter is to review theoretical foundations and some recent progress in this 
field. 

Most of the relevant physical phenomena are well known, such as interband tun¬ 
neling, Franz-Keldysh effect, Bloch oscillations, and Wannier-Stark localization. 
However, new experiments put these phenomena in a new context, which often leads 
to nontrivial observations, such as the generation of nonperturbative high-order har¬ 
monics in a solid due to Bragg-like scattering at the edges of the Brillouin zone 
[7, 8], a nearly instantaneous change in extreme-ultraviolet absorptivity and near- 
infrared refiectivity of a dielectric in the presence of a laser field as strong as several 
volts per angstrom [9], or the induction of electric curTent in an unbiased dielectric 
by similarly intense laser pulses [10]. Several decades of research on strong-field 
phenomena in solids and mesoscopic structures provide a solid ground for develop¬ 
ing new theoretical models adapted for new experimental conditions. At the same 
time, a description of extremely nonlinear processes that unfold during just a few 
femtoseconds requires approximations that may be different from those established 
for longer and less intense laser pulses. Both ab initio calculations that attempt to 
capture all the complexity of many-electron dynamics and “toy models” designed to 
deepen our understanding of basic phenomena are going to play an important role 
in extending our ability to control the optical and electric properties of solids with 
controlled light fields. 

The key parameters that determine the regime of strong-field light-matter inter¬ 
action with dielectrics and semiconductors are the Keldysh parameter 


COL\/mEg 
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( 1 . 1 ) 


the Bloch frequency 

eF]^a 
«B= , , 

and the Rabi frequency 

t 


( 1 . 2 ) 

(1.3) 


Here, TL is the amplitude of a linearly polarized electric field oscillating at an angu¬ 
lar frequency Ol, e > 0 is the absolute value of the electron charge, m is the reduced 
mass of an electron and a hole (m^* = a is the lattice period in the 

polarization direction of the electric field, and t/cv is the dipole matrix element re¬ 
sponsible for transitions between valence- and conduction-band states. Conditions 
that have to be fulfilled for perturbation theory to be applicable include )k ^ 1, 
Ob ^ Ol, ^ Ol, and «£g. Correspondingly, the main physical effects 
that make the interaction nonperturbative are interband tunneling (Yk ^1), Bloch 
oscillations (Ob > TZCOi), Rabi flopping (hQn > Eg), and carrier-wave Rabi flopping 
(f2R > Ol), the last two effects being particularly important for resonant excitations 
(fiOL « Eg). For a band gap of several electronvolts and a laser frequency in the 
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near-infrared spectral range (Al = 2ncl(0\^ ^ 1 pm), these effects become essential 
for > 1 V/A, which corresponds to a laser intensity of /l > 10*^ W/cm^ [11]. 


1.2 Main theoretical concepts 

For a theoretical description of phenomena that take place on a few-femtosecond 
time scale, it is common (although not necessarily correct) to neglect dephasing. In 
this case, the time-dependent Schrodinger equation (TDSE) 


ihd,\l/{t) =H{t)\i/{t) 


(1.4) 


fully describes nonrelativistic electron dynamics. In situations where phase and en¬ 
ergy relaxation processes are important, one has to use the more general formalisms 
of density matrices or nonequilibrium Green’s functions, which are beyond the 
scope of this chapter. In the following, we review some theoretical concepts de¬ 
veloped for the case of a single charged particle moving in a periodic potential. 
Furthermore, since the wavelength of visible or infrared light is much larger than 
the size of a unit cell, and electron velocities are much smaller than the speed of 
light, we can use the dipole approximation, which neglects the spatial dependence 
of the laser field while solving the TDSE: Fl = FL(f). 

The hrst step to solve the TDSE is to choose a gauge and a basis. The exact 
solution does not depend on this choice, but the chosen gauge and basis dictate ap¬ 
proximations that one may wish to make, and they influence the physical interpre¬ 
tation of results. In the dipole approximation, the two main options are the velocity 
and length gauges, which, in the following, will be abbreviated as “VG” and “LG”, 
respectively. These two gauges are related to each other by the following unitary 
transformation of the respective wave functions \j/\Q and i^lg^ 


V/vG(f)=exp --eAL(f)r V4-G(f) 
n 


(1.5) 


where 



( 1 . 6 ) 


is the vector potential of the laser field. 

The Hamilton operators in these two gauges take the following forms: 


Hlg = ^+t/(r)+eFL(f)r. 



(1.7) 


( 1 . 8 ) 


A big advantage of the velocity gauge for numerical simulations is that a homo¬ 
geneous external held does not destroy the spatial periodicity of the Hamiltonian, 
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SO that the Bloch theorem applies even in the presence of the field. At the same 
time, transformation (1.5) can be interpreted as a transition to a moving coordi¬ 
nate system, in which electrons acquire an additional momentum —eAL(f). In order 
to accurately account for such superhcial dynamics in a time-independent basis, a 
sufficient number of basis states and exact transition matrix elements are required. 
Calculations in the velocity gauge are particularly problematic in the limit of a static 
held (col —0). 

These problems are circumvented in the length gauge. The price for this is the 
fact that the term eF]^{t)r, which is responsible for the interaction with an exter¬ 
nal held, destroys the spatial periodicity of the Hamiltonian. Also, when periodic 
boundary conditions are applied in the length gauge, the interaction potential be¬ 
comes discontinuous at the boundaries of the unit cell. Nevertheless, once these 
difficulties are addressed, the length gauge becomes an appropriate choice for nu¬ 
merical simulations [12-16]. 

When the held of a laser pulse is nonresonant (/zCUl ^ Eg) and strong, it is fre¬ 
quently convenient to use a time-dependent basis of quantum states that adiabati- 
cally “adapt” themselves to the external held. In the length gauge, such an adiabatic 
basis is given by Wannier-Stark states. In the velocity gauge, this role is played by 
accelerated Bloch states, also known as Houston functions [17]. The following two 
subsections summarize the most important properties of these states. 


1.2.1 Wannier-Stark resonances 

For simplicity, let us consider a one-dimensional problem. For a constant external 
held parallel to the z-axis, the stationary Schrodinger equation reads 



The potential here is periodic; U{z + a) =U(z) with a lattice period a. Leaving the 
question of the existence of such eigenstates aside, one can ask which properties 
the solutions of (1.9) possess if they exist. From the periodicity of the potential, we 
immediately conclude that if (//(z) is an eigenstate with an energy E, then y/(z — a) 
is also an eigenstate with the energy E + eF^a. The additional term eFi^a is the 
energy required to move an electron against the laser held by one lattice period. This 
suggests that the eigenstates y/(z) should be localized functions. They are referred 
to as Wannier-Stark states. 

Wannier found approximate solutions to (1.9) by dehning an auxiliary problem 



( 1 . 10 ) 


where i is a band index, k is the crystal momentum, and Ei{k) is the energy of the 
unperturbed Bloch state. Wannier showed [18] that the solutions of this problem. 
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known as Wannier-Bloch states, decouple different bands. Specifically, knowing 
bj{z,k) that satisfies the periodic boundary condition in reciprocal space bi{z,k + 
2nja) = bi{z,k), an approximate solution of the time-dependent Schrodinger equa¬ 
tion with the Hamiltonian //lg can be written as 


Wiiz,t) =bi 



exp 





( 1 . 11 ) 


where the approximation consists in restricting the electron motion to a single band, 
that is, neglecting interband transitions. Obviously, |v/;(z,f)p is a periodic function 
of time, the period being equal to the period of Bloch oscillations: Tq = InjayR = 
iKtij{eFi^a). The Wannier-Stark states for a bulk crystal y/J^£^{z) are defined via the 
expansion 




( 1 . 12 ) 


A state Xj/Jg^ (z) is localized at a lattice site i. Explicit expressions for the Wannier- 
Stark states and their energies are 


= ^ ^ dkbi{z,k)e 

271 J-nfa 


-iiak 


=Ei + £eaFu 


(1.13) 

(1.14) 


where Ei is the mean energy of band i: 


Ei = 


a 



dkEi{k). 


The energies E^^^ form the so-called “Wannier-Stark ladder”—^plotted against 
TL, they are a set of straight lines, where the slope of each line dEj^^^/dFi^ = iea is 
determined by the lattice index £. 

The localization length of a Wannier-Stark state is given by 






(1.15) 


where A,- is the energy interval covered by band i. 

Our introduction to the Wannier-Stark states has so far followed the one given by 
Wannier [19]. The existence of Wannier-Stark states had been a subject of numer¬ 
ous disputes for three decades until they were experimentally observed in superlat¬ 
tices [20]. The core of these disputes was the question whether the Wannier-Stark 
states retain their physical significance when interband transitions are accounted 
for [21-23]. According to modern treatments [24-27], Wannier-Stark states should 
be viewed as resonances (metastable states) with lifetimes = l/f], and (1.14) 
should be generalized as 
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(1.16) 


Here, I] is the decay rate due to interband transitions, and 



(1.17) 


is called Zak’s phase [25, 28], where the integral is taken over a smooth closed path 
across the entire Brillouin zone, and ^,;(k) is the Berry connection or geometric 
vector potential [25-27]. An explicit expression for ^,y(k) is given by (1.26) in the 
next section. Zak’s phase plays an important role in the “modern theory of polariza¬ 
tion” [29]; it is equal to either 0 or n for crystals that possess inversion symmetry, 
and it can assume any value for other crystals. Direct measurements of Zak’s phase 
were performed for cold atoms in optical lattices [30]. 

Much of the mathematical complexity related to Wannier-Stark states is avoided 
in hnite systems, where the electron motion is restricted. In this case, it is common to 
refer to the exact length-gauge eigenstates of as Wannier-Stark states (without 
neglecting interband transitions by evaluating Wannier-Bloch states). These states 
have properties similar to those of the states introduced by Wannier. One of the most 
important differences is that the exact eigenstates of a Hamiltonian in one spatial 
dimension may not be degenerate. While the states dehned by (1.14) are strictly 
linear functions of Fl, and, for a certain value of Tl, some of these energies 
and may be equal to each other, the corresponding exact energies of a confined 
quantum system will have avoided crossings (anticrossings). An example of such a 
Wannier-Stark ladder is shown in Fig. 1.5. 


1.2.2 Accelerated Bloch states 

In section 1.2.1, we saw that the instantaneous eigenstates of the length-gauge 
Hamiltonian serve as a convenient basis for developing approximate solutions to 
the TDSE. In situations where the difficulties related to the length gauge outweigh 
its advantages, the velocity gauge may be a better choice for either numerical or 
analytical approximations, and the instantaneous eigenstates of the velocity-gauge 
Hamiltonian (1.7) may provide a more useful time-dependent basis. 

Let 0, k be a Bloch state with a band index i and a crystal momentum k: 



(1.18) 


In the coordinate representation, 0,,k(r) is a product of a plane wave and a lattice- 
periodic envelope function: 


</>!,k(r) = e''‘''M;,k(r) 


(1.19) 
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where (r + R) = Mi,k(r) for all R from the Bravais lattice. Let us now consider 
the instantaneous eigenstates of HyQ in the presence of a homogeneous external 
field: 

(^^^^^^ + U{r)^(p{t)=Eit)(pit). (1.20) 

Since the Hamiltonian is periodic in space, the Bloch theorem is applicable. Equa¬ 
tion (1.20) has the same form as (1.18), the momentum operator being substituted 
with p -l-eAL(f). The requirement that the solutions of (1.20) satisfy the Born-von 
Karman boundary conditions yields [31] 


<Pi.ko(r,0 =exp 




'/»r.k(r)(r), 


( 1 . 21 ) 


£,,ko(f)=£i(k(f)), 


( 1 . 22 ) 


where the time-dependent crystal momentum 

k(f)=ko+^AL(f) (1.23) 

n 

satisfies the acceleration theorem: hdk/dt = —eFi^{t). Here, ko is the initial crys¬ 
tal momentum, which the electron possessed prior to the interaction with the laser 
pulse. 

The states (r,0 are called accelerated Bloch states or Houston functions. We 
note that there is a close analogy between these functions and Volkov solutions of 
the TDSE [32], which are frequently used for describing the strong-held ionization 
of atoms and molecules [33]. 

To use Houston functions as a basis for solving the TDSE, it is convenient to 
make the ansatz 


i 

VAo(r,f) =L“''.ko(f)exp / df'£;(k(f')) 


^,-,ko(r,f), (1-24) 


which leads to the following system of differential equations [31, 34]: 

^a/.ko(0 = dt'AEij{k{t')) 

Here, A£’,y (k) = £’,(k) — ^^(k), and the matrix elements 

= (uk|iVk|7,k)^^ii = ^^dVM;k(r)VkM;,k(r) 


(1.25) 


(1.26) 


describe the optical transitions between bands, where the integration is performed 
over the volume of a unit cell. 

The relation between |,;(k) and the momentum matrix elements 
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p;j(k) = (/,k|p|7,k)^^ii = 


i/z 

a 


dVM*k(r)VrMj,k(r) 


(1.27) 


is given by 




moA£',j(k) ’ 


(1.28) 


where mo is the electron rest mass. This expression can only be applied as long 
as AEij{k) in the denominator is not equal to zero. The case of degenerate bands 
presents additional mathematical challenges [27, 35, 36]; in particular, the transition 
matrix elements ^,y(k) are singular at degeneracies [36]. 

For reference, we also give the relation between ^,j(k) and the matrix elements 
of the position operator between the Bloch functions [15, 37, 38]: 


(bk'|f|7,k)^= / dV([)*k/(r)r([)j.k(r)= i5,jVk + |;/k) 5(k-k'). (1.29) 


In the case a, ko (to) = Sij, where an electron is in band j before the external field 
is turned on, a simple approximate solution to (1.25) is [31] 







(1.30) 


provided that the excitation probabilities are small (|o!, ko |2 < 1 for i ^ j). 

This equation is a good starting point for numerous analytical approximations. 
In the case of a constant external field (Fl = const), it is convenient to rewrite the 
right-hand side as an integral over the crystal momentum k. Alternatively, the inte¬ 
gral over time can be approximately evaluated using the saddle-point method, which 
is especially useful in the case of a monochromatic external field. Such approxima¬ 
tions can be used to obtain analytical expressions for the rate of strong-field-induced 
transitions between valence and conduction bands of a dielectric or a semiconduc¬ 
tor. These transitions belong to the most important strong-held effects in solids, and 
they are discussed in the next subsection. 


1.2.3 Nonresonant interband transitions 

In 1928, Zener [39] used semiclassical arguments to show that a constant external 
held F makes valence-band electrons of a dielectric tunnel to the conduction band 
at a rate (per unit volume) 


J. _ e\F\a 

Tzener — /-» + ^Xp 

Inn 


2 eh\F\ 


The prefactor of the exponential function was found to be rather sensitive to a 
chosen method of approximation, but all such methods yield the same argument of 
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the exponential function. For example, Kane derived [40] the tunneling rate to be 
equal to 


.Hcane — ;; T/^ exp 


18 ;rr£', 


,2pl/2 


g 


2 eh\F\ ’ 


while Keldysh obtained [41] 


■Keldysh 


2|eF|V2,„i/4 

" 9;r2,372^5/4 


2 eh\F\ 


(1.31) 


It is common to state that such expressions for the tunneling rate have an “exponen¬ 
tial accuracy” [42]. 

If the external held is not constant but oscillating at a constant frequency, in¬ 
terband transitions may also occur as a result of absorbing a number of photons 
sufficient to overcome the band gap. It must be emphasized that there is no sharp 
distinction between nonresonant interband tunneling and multiphoton transitions. 
These are two asymptotic cases of interband excitations, which are distinguished by 
the Keldysh parameter (1.1); 'fk ^ 1 for multiphoton excitations and yK ^ 1 for 
tunneling. Equation (1.31) is valid in the latter case, where the laser held is strong 
and its frequency is small. In the intermediate regime (7 k ~ 1), it is impossible to 
unambiguously distinguish between contributions from multiphoton absorption and 
tunneling. This is reminiscent of the situation in atomic physics where the elec¬ 
tron motion under the potential barrier (tunneling) is known to be important even 
in the multiphoton regime [43]. Keldysh also derived a more general expression for 
the transition rate averaged over a laser cycle, which is applicable in all the three 
regimes [41]: 


fOKF = 


2®l 


OTOJl 


9k [ hp 


3/2 


2('m,A^)exp 


IJ 


E{a) 


(1.32) 


Q{n,N) = 


2K{a) 


1/2 


^exp 

«=0 


Kn[m)-Em 


X 


E{a) 

K^{2[N+1\-2N + N) 


2K{a)E{a) 


1 / 2 ' 


/ n -Eg ~ 2E{a) Eg 

Here, Eg is an effective ionization potential, the functions K{z) and E{z) are the 
complete elliptic integrals of hrst and second kind, ‘P{z) is the Dawson function, 
and [;icj denotes the integer part of x. We refer to (1.32) as the general Keldysh 
formula (GKF). 
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Keldysh parameter 


1.72 1.29 1.03 0.86 0.74 0.65 0.57 0.52 0.47 



Fig. 1.1 The probability that a 4-fs laser pulse excites a valence-band electron of Si02 to 
one of the conduction bands. The laser pulse, shown in the inset, has a central wavelength of 
Al = 800 nm. The dashed curves show probabilities evaluated by integrating analytical excitation 
rates: 1 — exp[— / F(f)dr]. For the curve marked as ’GKF’, the general Keldysh formula (1.32) was 
used as a cycle-averaged excitation rate; evaluating fGKF(0 we substituted Fl with the envelope 
of the laser pulse (dashed line in the inset). For the curve marked as ’tunneling’, we used the qua¬ 
sistatic approximation in the tunneling limit F(r) = rkeldysh(0 by substituting Fl with the electric 
field in (1.31). The solid curve is a numerical result obtained by solving the TDSE in a ID 
model [44]. The transition rates (1.31) and (1.32) were multiplied with constant factors to match 
the TDSE result for >k: < 1- 


Are these formulas, which were obtained decades ago for a monochromatic ex¬ 
ternal field, still useful in the case of few-cycle laser pulses? To address this ques¬ 
tion, we present, in Fig. 1.1, the outcomes of a simulation where a 4-fs 800-nm laser 
pulse interacts with a one-dimensional model medium that has properties resem¬ 
bling those of a-quartz [44]: a band gap of Fg = 9 eV, a lattice period of 5 A, and 
a reduced mass of m = 0.38mo. To compare the numerically evaluated excitation 
probabilities (solid curve) with the GKF predictions, we evaluated the excitation 
rate iGKF(0 using the real-valued pulse envelope in the place of the electric field 
Fl, which enters (1.32) via Y^^■ From rGKF(f)’ we estimated the excitation probabil¬ 
ity as fgkf 1 — exp[—/rGKF(f)d^] (thick dashed curve). The overall agreement 
is surprisingly good, given the fact that the laser pulse is shorter than two optical 
periods. Both the numerical and GKF results exhibit an oscillatory behavior, which 
appears due to closing and opening of multiphoton channels. It is analogous to the 
channel closing phenomenon in atomic physics [45-47]. The thin dashed curve in 
Fig. 1.1 represents the excitation probability evaluated with (1.31). Since Tkeidysh is 
the tunneling rate for a constant field, we used Fl)?) (rather than the pulse envelope) 
to evaluate the excitation probability by the laser pulse. This procedure is known as 
the quasistatic approximation [48]. In this example, the tunneling formula (1.31) is 
inaccurate for Fl < 1.2 V/A, which corresponds to tk 0.9. 

The pioneering work by Keldysh was followed by numerous investigations. A 
few that we would like to point out here are analytical results obtained using the adi¬ 
abatic approach, where parity selection rules were accounted for [34,49], derivation 
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of expressions for arbitrary A^-photon transition probabilities [50], and development 
of the Keldysh-like theory for cosine-shaped bands [51]. 


1.3 Strong-field-driven electron dynamics in crystals 
1.3.1 A numerical example 

A strong electric held drives interband transitions, it accelerates charge carriers, and 
it can also cause transient changes in the optical properties of a medium without 
necessarily exciting electrons to conduction bands. Combined with controlled opti¬ 
cal helds, these three classes of physical phenomena enable the manipulation of the 
electric and optical properties of a medium over time intervals much shorter than a 
period of optical oscillations. The main examples of such controlled helds are laser 
pulses with the stabilized carrier-envelope phase (CEP) and optical waveforms [52]. 

We begin this section by illustrating such effects in a simulation where the TDSE 
was solved in one spatial dimension for a periodic potential. In Eig. 1.2, we show 
obtained by solving the Houston-basis equations (1.25) and representing 
the conduction-band populations in the extended-zone scheme, where the range of 
crystal momenta covered by the n-th conduction band (n > 1) is n — 1 < |k| /kmax < n 
with kmax = Ttla. According to the acceleration theorem, the ballistic motion of an 
electron wave packet is described by hd{k)/dt = —eFi^it), where {k) is the mean 
crystal momentum. A continuous change of {k) at the borders |k|/kmax = « corre¬ 
sponds to transitions between different conduction bands. Such interband transitions 
are closely related to Landau-Zener transitions, as they occur at crystal momenta 
where the energy gap between two adjacent bands is particularly small. If a charge 
carrier remains within its current band as it crosses a Brillouin zone edge, its crystal 
momentum changes abruptly. This can be interpreted as a Bragg-like reflection of 
an electron wave off the crystal lattice. Eig. 1.2 shows that both Bragg reflections 
and interband transitions play an important role when the electron motion is driven 
by a near-infrared held. 

In Fig. 1 .2(a), the electron wave packet is initially placed in the lowest conduction 
band of the ID model of a solid, and even though the laser pulse with —1.6 pm is 
strong enough to accelerate electrons out of the first Brillouin zone, its amplitude is 
insufficient to induce transitions between conduction and valence bands. The figure 
illustrates the importance of transitions between different conduction bands, and 
shows that they predominantly occur atk — nn/a, n S Z. 

For the simulations presented in Fig. 1 .2(b), all electrons were initially placed in 
the valence bands. The laser pulse had the same peak value of the vector potential 
as in Fig. 1 .2(a), but it had a shorter wavelength Al = 800 nm, so that the electric 
held of the pulse was twice as strong, and it was strong enough to excite electrons 
from the uppermost valence band. One can see that the population of conduction 
bands does not constantly increase with time, as one would expect from rate models 
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Fig. 1.2 The time-dependent population distribution of conduction-band electrons in the pres¬ 
ence of a few-cycle laser pulse. These plots were obtained by solving (1.25) for a ID model 
of Si02 with the same lattice potential as that used in [44]. The normalized conduction-band 
population {maxn(t,k) = 1) is plotted against time and crystal momentum k in the extended- 
zone scheme, where the range of crystal momenta covered by the n-th lowest conduction band 
is fi — 1 < |k|/^,„ax < « with kmax = Tt/a. The dashed curves show AL/A:niax- (a) An electron wave 
packet is initially placed in the lowest conduction band (o£„uj(-(fo) = exp[—an 8-fs 
laser pulse with a central wavelength of Al = 1.6 pm and a peak field of Fl = 0.7 V/A is sufficiently 
strong to accelerate the electrons out of the first Brillouin zone, (b) All electrons are initially in the 
valence bands; a 4-fs pulse with Al = 800 nm and Fl = 1.4 V/A excites them to the conduction 
bands, where they are accelerated by the laser field. 


discussed in section 1.2.3. Instead, there is a transient increase of the population at 
the extrema of the electric field (at the zero crossings of the vector potential). Such 
excitations are called “virtual” [53] because they only exist as long as the exter¬ 
nal field is present. If a laser pulse is so weak that the final excitation probability 
can be neglected, such virtual excitations represent the distortion of bound states 
of the crystal—the distorted valence-band states have nonzero projections onto the 
field-free conduction-band Bloch functions. In the strong-held regime, no clear dis¬ 
tinction between the “virtual” and “real” excitation can be made. Nevertheless, in 
both weak- and strong-held regimes, the transient increase of conduction-band pop¬ 
ulation has an effect on experimentally observable quantities [53, 54]. 

Currently, there is no measurement technique that would reveal all the details 
of strong-held-driven electron dynamics like those illustrated by Fig. 1.2. However, 
time-resolved measurement techniques do provide indirect access to this informa¬ 
tion. Some evidence of subcycle dynamics in the strong-held excitation of electron 
in Si02 was presented in [55] by measuring the polarization rotation of an ellipti- 
cally polarized pulse transmitted through a thin glass plate, as well as in [56] using a 
noncollinear pump-probe measurement scheme. Bragg-like scattering of electrons 
was found to contribute to the generation of nonperturbative high-order harmonics 
in solid samples [7, 8]. The reversible held-induced change of absorption in the ex¬ 
treme ultraviolet spectral range was observed by probing the effect of an intense 
near-infrared held on a thin silica plate using an attosecond pulse of extreme ultra- 
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violet radiation as a probe [9]. A subcycle turn-on of electric current in a dielectric 
and its manipulation with CEP-stabilized pulses was demonstrated by measuring the 
residual polarization induced by the laser light [10]. This last effect is most closely 
related to the topic of this chapter, so it is discussed in more detail in the following 
section. 


1.3.2 Ultrafast injection and control of current in dielectrics 

According to Fig. 1.2(b), an optical field that is strong enough to excite valence-band 
electrons of a dielectric to its conduction bands is also strong enough to significantly 
accelerate the created charge carriers, thus driving electric current. One of the most 
important hndings in [10] was that such electric current can be turned on within 
a fraction of a half-cycle of a short intense laser pulse. This was demonstrated by 
irradiating a Si02 sample placed between two gold electrodes with CEP-stabilized 
laser pulses (see Fig. 1.3) and measuring the current induced by the pulses in an 
external circuit. It was found that a short laser pulse with the electric field directed 
perpendicularly to the electrodes was able to leave the sample in a polarized state, 
implying that a certain electric charge was displaced by the pulse. By varying the 
CEP of the pulse, it was possible to control the amount of the displaced charge. 
This fact alone already suggests that the observed effect should be controlled by the 
electric field (rather than the envelope) of the laser pulse, but the most convincing 
evidence for the subcycle turn-on of electric current was provided by pump-probe 
measurements, where the sample was irradiated by a pair of pulses; an intense “in¬ 
jection pulse” polarized parallel to the electrodes and a relatively weak “drive” pulse 
polarized perpendicularly to them. By observing how the displaced charge depends 
on the delay between the two pulses, it was possible to conclude that the injection 
pulse makes the Si02 sample conductive within a time interval < 1 fs. These mea¬ 
surements were well reproduced by simulations: the one-dimensional tight-binding 
simulations in the original paper [10], a model that used a one-dimensional pseu¬ 
dopotential [44], and recent ab initio three-dimensional simulations [57]. Neverthe¬ 
less, these observations permit several interpretations, which we present in the rest 
of this section. 


Semidassical interpretation 

The most intuitive interpretation is based on the observation that the strong-held- 
driven motion of an electron wave packet in the conduction bands of a dielectric 
largely obeys the acceleration theorem (see Fig. 1 .2(a))— apart from interband tran¬ 
sitions and occasional Bragg scattering, which splits and reunites electron wave 
packets, each wave packet moves as a classical particle with an effective mass that 
depends on the mean crystal momentum. This suggests that the residual polariza¬ 
tion induced by the laser pulse may be interpreted in semidassical terms. A rigorous 
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Fig. 1.3 A schematic representation of the experimental arrangement used in [10], An intense 
few-cycle laser pulse induces electric current in the dielectric (SiOa) placed between two gold 
electrodes. The measured signal is proportional to the net charge displaced by the laser pulse, and 
it is conh'olled by the carrier-envelope phase of the laser pulse. 


approach to this interpretation consists in writing the current density averaged over 
a unit cell 

in the basis of Houston functions, where the contribution from an electron with an 
initial crystal momentum kp is given by 


jko (0 = -eY,\ (0 1^ Vi (k(f)) 

i 

2e f 

-E a-k„ (?)«;■ ko(OPii(k(f)) exp 

ij<i f 



(1.34) 


Here, k(f) is dehned by (1.23), mo is the electron rest mass, P!y(k) are the mo¬ 
mentum matrix elements (1.27), and v,(k) = Vk£’,(k)//! = p, 7 (k)/mo is the group 
velocity in band i. The hrst sum on the right-hand side of (1.34) is responsible for 
the current due to the ballistic motion of charge carriers. In the semiclassical in¬ 
terpretation, the contribution from this term to the residual polarization (displaced 
charge density) 

/ ^max 

J(f)d? 

is assumed to be much larger than that from the second sum, which describes effects 
related to interband coherences. 

In the tunneling regime (Yk. ^1), charge carriers are predominantly created at 
the extrema of the electric held, each of which launches an electron wave packet. 
For a wave packet launched at a time fp, it is convenient to introduce a semiclassical 
displacement. 


s{to)= f v(k(f))df 
Jto 


(1.35) 
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where the group velocity v(k) should correspond to the most probable quantum 
path of the wave packet in reciprocal space. As long as k(f) is not limited to the first 
Brillouin zone, this approach is most useful if the probabilities of Bragg scattering 
at the edges of the Brillouin zone are either negligibly small or close to 100%. The 
contribution from each wave packet to the final polarization P(fmax) is the product 
of the charge carried by the wave packet and its semiclassical displacement s(fo). 

This kind of semiclassical analysis explains some outcomes of numerical simu¬ 
lations [58]. In particular, it explains the observation that, for moderate laser inten¬ 
sities, the residual polarization scales as P p 2 N+\ where N — E^/{hcO\_) is 
the ratio of the band gap to the photon energy. As long as the Keldysh parameter is 
sufficiently large to view interband excitations as the result of absorbing N photons, 
the excitation probability scales as p being the peak laser intensity, 

while the semiclassical displacement is proportional to Pl- Thus, the product of the 
charge and the displacement is proportional to 


Interference of multiphoton pathways 

For moderate laser intensities, where perturbation theory is expected to yield at least 
qualitatively correct predictions, it is also possible to interpret the optically con¬ 
trolled electric current in terms of interference between different multiphoton exci¬ 
tation pathways. This interpretation is a generalization of ideas developed in the field 
of coherent control, where irradiating a semiconductor by two monochromatic laser 
beams with frequencies coi and 0)2 = 2coi was found to induce an electric current 
sensitive to the relative phase between the two beams [59]. This phase sensitivity is 
due to the interference between single- and two-photon absorption processes. The 
density n(k) of electrons excited to the conduction band at a crystal momentum k is 
determined by the phase parameter = 2(po,j — (pay,. Furthermore, n(k) is, in gen¬ 
eral, an asymmetric function of k because the transition amplitudes for the one- and 
two-photon channels have different symmetries with respect to the transformation 
k —) —k [59]. Injecting currents through interfering photoexcitation pathways was 
investigated in theory and experiments for semiconductors [59-66] and molecular 
wires [67]. 

In the case where charge carriers are excited by a few-cycle pulse with a central 
frequency that is much smaller than the band gap, the interfering photoexcitation 
pathways are multiphoton excitation channels. It is convenient to analyze them in 
the velocity gauge using the basis of Bloch states, where an external homogeneous 
electric field only induces transitions between states with the same the crystal mo¬ 
mentum. Due to the large bandwidth of ultrashort laser pulses, it is possible to make 
the same transition by absorbing different numbers of photons, as it is schematically 
shown in Fig. 1.4. 

Within this picture, the scaling law Q oc ^ 2^+1 interpreted as a result of 

interference between quantum pathways that involve absorbing N and N +1 pho¬ 
tons. The corresponding probability amplitudes are proportional to and 
respectively. When the two pathways interfere, the excitation probability, which is 
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Fig. 1.4 Current injection via interference of multiphoton excitation channels. Charge carriers can 
be injected from a valence band (VB) into a conduction band (CB) via different multiphoton excita¬ 
tion pathways if the spectral bandwidth of the applied laser light is sufficiently large. For example, 
the excitation at crystal momenta k 2 and —k 2 may be the outcome of absorbing either three high- 
energy or four low-energy photons from the same ultrashort laser pulse. Quantum interference 
between the odd- and even-numbered contributions in each pathway determines the population in 
the final state. An asymmetry in the conduction-band population {n(k) ^ n(—k)) results in a net 
current density is formed inside the material. 


the squared modulus of the sum of the probability amplitudes, is described by an 
expression that contains the product of the two amplitudes. This product, which is 
proportional to determines the induced electric current. Even though these 

arguments by no means form a rigorous proof, and these considerations are only 
applicable in the multiphoton regime, the interpretation in terms of interfering mul¬ 
tiphoton channels was shown to be a plausible one, explaining not only the scaling 
with intensity, but also the fact that a CEP-sensitive displaced charge can only be 
observed if the laser pulse is sufficiently broadband [44]. 


Adiabatic metallization 

Another interpretation of the optical-held-induced current in dielectrics accompa¬ 
nied the hrst publication of these experimental results [10]. This mechanism re¬ 
lies on the effect of “adiabatic metallization”, predicted theoretically for dielectric 
nanohlms [68, 69]. It was found that a strong electric held can signihcantly and 
reversibly change the optical and electric properties of a sufficiently thin dielectric 
nanohlm. During the interaction with the held, the real part of the dielectric con¬ 
stant may even become negative, which is a property attributed to metals. These 
effects are best understood in the basis of instantaneous eigenstates of the length- 
gauge Hamiltonian, which are the Wannier-Stark states discussed at the end of 
section 1.2.1. Even though these results were obtained for nanohlms, similar ef¬ 
fects may be expected in bulk solids [70], provided that the localization length of 
Wannier-Stark states does not exceed a few lattice sites, which is indeed the case 
for held strengths EL > 1 V/A. 
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Electric field (V/A) 


Fig. 1.5 Eigenstates of the Wannier-Stark Hamiltonian (1.8) for a 50-nm nanofilm of Si02 in 
dependence of the electric field. The labels Al = i and Al = 2 indicate the anticrossings that 
correspond to Wannier-Stark states being localized one and two lattice sites apart, respectively. 
Closed circles represent occupied states, while open circles represent unoccupied states. The figure 
is adapted from Ref. [71]. 


Using the Wannier-Stark states as a time-dependent basis, the interaction with 
an intense pulse can be analyzed in terms of adiabatic and diabatic transitions at 
avoided crossings. This analysis [10, 71] shows that the strong field may cause 
a reversible quantum transition from an insulating state to a state with an in¬ 
creased conductivity. This transition occurs when the energy gaps at avoided cross¬ 
ings between Wannier-Stark states become sufficiently large for electrons to adi- 
abatically pass such an anticrossing (see Fig. 1.5). The anticrossing gap takes its 
largest value when the relevant Wannier-Stark states are localized at neighboring 
lattice sites, which happens when the electric field in the medium is as strong as 
TcUt — A'g/This condition is equivalent to /zCUn — Ag. For Si02, ~ 3 V/A, 

which is above the damage threshold, but even the penultimate anticrossing, which 
occurs at Ferit « 1.5 V/A, was found to have a sufficiently large energy gap. This in¬ 
terpretation of optical-field-controlled current as a result of adiabatic metallization 
was also supported by the observation of the optical-field-induced transient reflec¬ 
tivity in the strong-held regime [9]. 


1.4 Summary and outlook 

The main message of this chapter is that the control over the electric held of light 
pulses can be transformed into the control over strong-held-driven electron dynam¬ 
ics, which may open attractive opportunities for both basic research and applications 
in ultrafast signal processing. One example of such an application is a solid-state 
CEP detector [72], and more advanced applications may emerge in the future. At- 
tosecond science provides powerful tools and techniques for studying strong-held 
dynamics in solids, which may lead to important insights into fundamentally impor- 
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tant phenomena. So far, only the first few steps have been done in this direction, and 
many questions remain open. Further research is necessary to verify conjectures 
made in the interpretation of recent measurements, and it is yet to be understood 
how electron-electron interaction, scattering, and dephasing phenomena affect our 
ability to launch and manipulate electron wave packets. One of the most important 
goals is to identify those extremely nonlinear effect that are largely reversible on a 
few-femtosecond time scale, as only such effect may serve as a basis for applica¬ 
tions in signal processing. On the experimental side, there is a large potential for 
exploiting various laser sources, intense optical waveforms, and (nano) structures 
designed for strong-held measurements. 
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